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MHD  PLASMA  CHANNEL  RESPONSE  TO  PROPAGATING  ION  BEAMS 
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Mega-ampere  beams  of  1-2  MeV  protons  and  deuterons  have  recently  been 
2 

extracted  from  100  cm  Pinch  Reflex  diodes  coupled  to  high-power  trans- 

1  2  12 
mission  line  generators.  Focusing  to  .4  MA/cm  has  been  reported.  ’ 

These  light  ion  beams  are  appropriate  for  inertial  confinement  fusion 

experiments  when  a  large  number  of  them  are  individually  transported 

3  2 

several  meters  and  combined  onto  a  pellet.  To  this  end,  100  kA/cm  beams 
have  been  transported  a  meter  distance  through  50  kA  discharges  in  1  torr 

4 

air.  The  azimuthal  magnetic  field  produced  by  the  discharge  current 
confines  the  charge-  and  current-neutralized  beam  within  the  discharge 
diameter.  In  the  present  theoretical  work,  the  MHD  response  of  transport 
channels  to  ignition-system-level  beams  is  investigated.  Simple  analytic 
estimates  for  channel  response  agree  with  a  1-D  radial,  2  temperature  MHD 
code  for  moderate  density  changes.  Results  indicate  that  for  3  to  5  MeV 
protons,  mega-ampere-level  beams  can  be  transported. 

The  discharge  current  1^  required  to  transport  a  beam  with  given 
injection  conditions  is  determined  from  conservation  of  pz~conical  momentum. 
For  an  ion  entering  a  discharge  at  radius  rQ  with  speed  vq  at  angle  0  to  the 
z  axis  with  no  angular  momentum, 

vz(r)  *  vQ  cos  0  +  qi/“i Jr  B(r)dr  (1) 

r 

o 
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where  B  is  the  azimuthal  magnetic  field.  For  proton  beams  and  B(r)  ^  r, 

I  u(A)  a  .5  10  J  8  2  v  (cm/s)  where  9  (rad)  is  the  maximum  injection  angle 
cn  mo  m 

For  50  kA  discharge  currents  and  several-MeV  protons,  this  expression 
requires  9m  ^  10°.  A  uniform  current  distribution  is  appropriate  for  few- 
electron-volt,  centimeter-diameter  discharges  established  microseconds 
before  beam  injection.  Once  the  beam  is  injected  into  the  discharge, 
kinetic  pressure  and  magnetic  forces  on  the  plasma  return  current  induced 
by  beam  injection  drive  channel  expansion  on  the  50  ns  beam  duration  time 
scale.  The  expanding  plasma  convects  magnetic  field  outward,  thereby 
increasing  the  beam  radius.  At  the  same  time,  the  beam  experiences  energy 
losses  due  to  collisional  deposition  and  deceleration  in  induced  electric 
fields.  The  plasma  channel  density  must  be  low  enough  to  prevent  excessive 
energy  losses  but  high  enough  to  prevent  excessive  expansion. 

In  order  to  study  the  response  of  the  channel  to  injected  beams  with 
a  1-D  radial  code,  azimuthal  symmetry  and  small  axial  gradients  are  assumed 
The  equations  solved  in  the  code  include  the  MHD  equations  with  classical 
electron  transport  coefficients  and  the  appropriate  Maxwell's  equations. 
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In  the  above,  j  and  j^  are  the  plasma  and  beam  current  density,  is  the 
electron-ion  equipartition  time  and  Prad  is  the  optically-thin  radiation  loss 
rate. 
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The  standard  channel  prior  to  beam  injection  is  6  x  10  g/cm 

deuterium  (n.  ■  2  *  10^8  cm  "*) ,  with  I  ,  =  50  kA,  T  =  T .  w  3  eV  inside  a 
l  cn  e  i 

channel  of  radius  r£h  *  0.56  cm.  Deuterium  is  chosen  for  simplicity  of 
atomic  physics  and  radiation  processes,  and  for  collisional  beam  losses  lower 
than  a  comparable  mass  density  of  hydrogen.  The  plasma  channel  is  surrounded 
by  a  cold  gas  blanket.  Radial  profiles  for  the  channel  prior  to  beam 
injection  (Fig.  1)  are  derived  from  pressure  balance  for  1^  uniformly 
distributed.  The  transported  beam  is  modeled  by  j^  *  j^^t)  ^or  r  ^  ^rch’ 
and  Jb  “  jbo  (1  “  r/  rcli)/(l  -  y)  for  yrch  <  r  4  rch  with  .1  ><  y  .<  .35.  For 
beam  duration  t,  j  -  j^  (1  +  t/r)/2.  Ion  energyC^  increases  in  time 


3 
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like  j.,  reaching Cbm  at  t  *  x.  This  energy  ramp  permits  beam  power 
multiplication  by  axial  bunching  during  transport. 

In  Eq.  3,  the  magnetic  expansion  forces  usually  dominate  over 
kinetic  pressure.  Assuming  constant  radial  acceleration  due  to  j^B/c, 
plasma  expansion  is  approximated  by 

X  i  3  *  10  2  (A/cm2)  jch  (A/cm2)  t2 (s) /p (g/cm3)  (13) 


where  a  fluid  element  initially  at  radius  r  is  displaced  to  (1  +  x)r  at 

time  t,  i  ,  is  the  discharge  current  density  and  j.  -j  is  assumed  for 
cn  b  p 

jcb  <<  j^.  Radial  channel  and  beam  expansion  will  be  small  when  x  <<  1. 

2 

The  standard  channel  parameters  give  x  —  1/2  after  50  ns  for  a  1  MA/cm 
beam. 


Ion  compressional  heating  in  Eq.  4  is  always  less  important  than  the 

electron-heating  mechanisms  on  the  right-hand  side  except  at  the  channel 

2 

boundary.  For  =  5  MeV  and  =  1  MA/cm  ,  collisional  deposition  dominates 

over  return-current  heating  above  3  eV.  The  proton  stopping  power  in 

3  2 

deuterium  is  then  S  =  10  /£  MeV/g/cm  .  The  associated  temperature  rise 

b 

3 

in  the  plasma  is  6(T  +  T.)  'v  7  x  10  i, t/C,  and  is  of  order  100  eV  for 

e  i  b  b 

the  above  beam  after  50  ns.  For  2.3-meter  transport,  collisional 

beam  energy  losses  are  less  than  20%  of  the  initial  beam  energy  when 

P  <  9  x  lo  7G  2.  This  criterion  combined  with  that  for  inertial  confinement 
b 

(Eq.  13  with  x  <  1/2)  yields  7  *  10  S  P  ^  9  "  10  ^C^2,  i,e*  ^"b  >  2‘®  MeV* 
Energy  losses  due  to  deceleration  in  the  axial  electric  field  of  Eq.  8  must 
also  be  considered.  Early  in  time,  E  is  established  by  return-current  flow 
in  the  resistive  plasma.  Later,  high  plasma  conductivity  and  expansion 
velocity  cause  the  hydrodynamic  term  to  dominate.  Energy  losses  due  to 


« 


late-time  electric  field  are  comparable  to  collisional  losses  for  the  above 
density  regime. 

2 

Figure  2  shows  the  channel  at  50  n&  for  j  =  .6  MA/cm  ,  and 

18 

C  =  5  MeV.  Expansion  has  reduced  the  initial  density  to  1.2  -  1.4  *  10 

_3 

cm  .  This  agrees  with  the  calculated  value  of  x  =  .3  from  Eq.  13.  A 

figure  of  merit  for  ion  confinement  is  defined  by  the  integral  in  Eq.  1 

extended  from  0  to  r  ,  .  The  code  calculates  a  reduction  in  its  value  from 

ch 

an  initial  5  kG-cm  to  3*3  kG.cm  for  Fig.  2  in  agreement  with  simple  scaling. 

At  the  conditions  shown,  transported  ions  loose  about  .14  MeV/m  collision- 

ally  and  about  .1  MeV/m  due  to  electric-field  deceleration.  Early  in  time, 

when C.  =  2.5  MeV,  corresponding  values  are  .28  MeV/m  and  .03  MeV/m. 
b 

Figure  3  illustrates  the  response  of  the  standard  channel  to  a  beam 

2 

of  comparable  power  density  as  in  Fig.  2  but  with  =  2.75  MA/cm  and 
Cbm  =  3  MeV.  This  case  is  very  different  from  the  preceding  one  since  the 
channel  has  had  time  to  expand  and  to  bounce  from  the  cold  gas  blanket.  The 
structure  of  the  various  fluid  quantities  is  much  more  complicated  and 
indicates  that  different  physics  applies.  For  example,  instabilities  might 
play  an  important  role  which  would  change  the  results  quantitatively. 
Estimates  from  Eq.  13  can  still  be  made,  but  only  within  a  radius  of  about 
2  mm  where  plasma  expansion  still  occurs.  In  that  region,  the  code  is  in 
agreement  with  simple  estimates.  An  interesting  feature  of  these  results 
is  that  the  vB/c  term  gives  an  electric  field  which  will  accelerate  the  ion 
beam  if  it  exceeds  the  collisional  slowing  down  (ohmic  electric  fields  are 
small  due  to  high  conductivity). 

In  addition  to  the  case  of  the  unbunched  beam  which  prevails  at  the 
beginning  of  the  channel,  the  plasma  channel  response  to  a  bunched  beam  near 
the  exit  of  the  channel  has  been  investigated.  It  is  found  that  a  beam  of  a 
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given  current  and  voltage  is  more  easily  transported  in  its  bunched  state  due 
to  the  reduced  hydrodynamic  expansion.  This  can  be  seen  from  Eq.  13  with 
jfet  constant  so  that  x  depends  only  on  For  the  10  ns  pulse  duration 

of  a  bunched  beam,  the  radial  expansion  is  5  times  smaller.  Alternatively, 
the  density  can  be  chosen  smaller  to  maintain  a  constant  x*  For  constant 
collisional  losses  in  transit  through  the  channel  (i.e.  /pdz  =  const.),  the 
density  at  the  beginning  of  the  channel  can  then  be  increased  to  reduce  the 
channel  expansion  there. 

In  conclusion,  the  hydrodynamic  plasma  response  to  light-ion  beams 
has  been  studied  using  simple  estimates  and  numerical  solutions.  For  3-5  MeV 

proton  beam  propagation,  deuterium  z-discharges  were  studied  for  injected 

2  —5  3 

MA/cm  ,  50  ns  beams.  Discharge  channel  densities  'v  10  g/cm  were  found 

to  permit  transport  of  these  beams  in  both  unbunched  and  bunched  states  over 

the  several-meter  distances  required  for  ICF  ignition  experiments. 

This  work  was  supported  by  the  Department  of  Energy. 
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